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Abstract 

Background: Pseudorabies virus (PRV), an alplia-lierpesvirus of swine, is a widely used model organisnn in 
investigations of the molecular pathomechanisms of the herpesviruses. This work is the continuation of our earlier 
studies, in which we investigated the effect of the abrogation of gene function on the viral transcriptome by 
knocking out PRV genes playing roles in the coordination of global gene expression of the virus. In this study, we 
deleted the usl gene encoding the ICP22, an important viral regulatory protein, and analyzed the changes in the 
expression of other PRV genes. 

Results: A multi-timepoint real-time RT-PCR technique was applied to evaluate the impact of deletion of the PRV 
usl gene on the overall transcription kinetics of viral genes. The mutation proved to exert a differential effect on 
the distinct kinetic classes of PRV genes at the various stages of lytic infection. In the usl gene-deleted virus, all the 
kinetic classes of the genes were significantly down-regulated in the first hour of infection. After 2 to 6 h of 
infection, the late genes were severely suppressed, whereas the early genes were unaffected. In the late stage of 
infection, the early genes were selectively up-regulated. In the mutant virus, the transcription of the iel80 gene, the 
major coordinator of PRV gene expression, correlated closely with the transcription of other viral genes, a situation 
which was not found in the wild-type {wt) virus. A 4-h delay was observed in the commencement of DNA 
replication in the mutant virus as compared with the wt virus. The rate of transcription from a gene normalized to 
the relative copy number of the viral genome was observed to decline drastically following the initiation of DNA 
replication in both the wt and mutant backgrounds. Finally, the switch between the expressions of the early and 
late genes was demonstrated not to be controlled by DNA replication, as is widely believed, since the switch 
preceded the DNA replication. 

Conclusions: Our results show a strong dependence of PRV gene expression on the presence of functional usl 
gene. ICP22 is shown to exert a differential effect on the distinct kinetic classes of PRV genes and to disrupt the 
close correlation between the transcription kinetics of iel80 and other PRV transcripts. Furthermore, DNA replication 
exerts a severe constraint on the viral transcription. 
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Background 

The pseudorabies virus (PRV), an alpha-herpesvirus, is 
the etiological cause of Aujeszkys disease of swine [1]. 
PRV is related to the human pathogen varicella-zoster 
virus (VZV) and herpes simplex virus types 1 and 2 
(HSV-1 and -2), and the animal herpesvirus bovine her- 
pesvirus type 1 (BHV-1). PRV is widely used as a model 
organism in investigations of the molecular patho- 
mechanisms of the herpesviruses [2], and is a useful tool 
for the mapping of neural circuits [3,4]. Attempts have 
additionally been made to utilize this virus as a gene de- 
livery vector [5,6] and an oncolytic agent [7]. Besides the 
lytic phase, alpha-herpesviruses can enter a latent state, 
where they transcribe a limited set of c/5-antisense RNAs 
[8]. Traditionally, the lytic herpesvirus genes are classi- 
fied into three kinetic categories: immediate-early (IE) 
genes, early (E) genes and late (L) genes. On a finer 
scale, an intermediate category, the early/late (E/L = 
delayed early) genes can also be distinguished [9]. PRV 
encodes a single IE gene, the iel80 gene. The IE 180 pro- 
tein, a transactivator, is the principal coordinator of the 
overall gene expression of the virus. E genes encode 
proteins required for the nucleotide metabolism and 
DNA replication. Other E genes such, as the early pro- 
tein 0 {epO) [10] and ul54 genes [11], encode transcrip- 
tional regulators. Most of the L genes code for structural 
elements of the virus. ICP22 is one of the five IE 
proteins of HSV-1, which is encoded by the usl gene. 
Intriguingly, a large part of the HSV usl gene is located 
in the unique US region, whereas its promoter and a 
short 5' portion of the transcribed region are in the 
inverted repeat (IR) segment. In PRV, however, the entire 
usl gene (earlier called the rsp40 gene) resides in the IR 
region; this gene is therefore represented in two copies 
in the PRV genome. There is no consensus as to whether 
the PRV usl gene is expressed in IE [12] or E kinetics 
[13]. We demonstrated in an earlier analysis that this 
gene is expressed in atypical kinetics, and that it is obvi- 
ously not an IE gene [9]. The function of the ICP22 
polypeptide has primarily been analyzed in HSV-1. The 
investigations have revealed that ICP22 is a multifunc- 
tional protein that plays roles in various aspects of HSV 
pathogenesis. It has not yet clearly established whether 
ICP22 acts to repress E genes [14] or to enhance the 
transcription of L genes [15]. It has been shown that 
not all L genes require this transactivator for their 
expressions [16]. The BICP22 protein of BHV-1, a 
homologue of ICP22, has been demonstrated to exert a 
general repressive effect on each kinetic class [17]. Rice 
and coworkers reported that ICP22 acts at the level of 
transcriptional regulation [18]. However, the level of 
ICPO mRNA was also reduced in the usl knockout (KO) 
HSV [16], which raises the question of whether the dir- 
ect cause of the reduced transcription is the lack of usl 



gene activity or the low ICPO mRNA level. ICP22 has to 
be phosphorylated by the viral UL13 protein kinase in 
order to accomplish the transcriptional activation of L 
genes [19]. An additional function of the ICP22 polypep- 
tide is associated with the alteration of the activity of 
cyclin-dependent kinase cdc2, a regulator of the cell cycle, 
which results in a selective up-regulation of HSV L genes 
during lytic infection [20]. Furthermore, HSV ICP22 also 
acts to modify the phosphorylation of RNA polymerase II 
(RNAP II) [21], which carries out the transcription of viral 
genes. One of the major control regions of RNAP II is its 
carboxy terminal domain (CTD), residing on the large 
subunit of the molecule. The CTD, containing multiple 
repeats of a heptapeptide sequence, serves as a binding 
site for various cellular proteins involved in the regulation 
of transcription. ICP22 is presumed to trigger the loss of 
Ser-2 phosphorylation on the CTD, and thereby modif)^ 
the activity of RNAP II [22]. A novel function of ICP22 
was recently identified, involving alteration of the chap- 
eron localization of the host cells [23]. It has been shown 
that ORF63, the ICP22 homolog of VZV, does not alter 
RNAP II phosphorylation and the host chaperon machin- 
ery [22], which might indicate that ICP22 acts in a 
species- or genus-specific manner. In the present study, 
we have investigated the effects of usl gene deletion on 
the overall transcription of PRV genes. 

Results and discussion 

Experimental design 

An insertion mutant PRV strain was constructed which 
contains the mutation in both copies of the usl gene. 




time (h) 

— •— low titer (MOI: 0.1 pfu/cell) infection 
■ -O' ■ high titer (MOI: 10 pfu/cell) infection 

Figure 1 Transcription Icinetics of the usl gene of the wt virus 
in low-MOl (0.1 pfu/ceil) and high-MOl (10 pfu/cell) infection 
conditions. The low-MOl experiment was conducted only up to 6 h 
pi: after this time, the newly generated virions may be released from 
the infected cells and initiate new infections in the originally non- 
infected cells. The Figure shows a 2-h delay in the initiation of the 
steep rise in gene expression. 
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Mutation of a DNA sequence in the internal repeat re- 
gion is copied to the terminal repeat by a mechanism 
called equalization [24]. From among the 70 PRV genes, 
32 were selected for the transcription analysis, which 
reside at the upstream position of the tandem gene 
clusters and which represent each kinetic class of PRV 
genes. The reason for this choice was to exclude the 
distorting effect of the transcriptional read-through 
exerted by the upstream genes on the downstream 
genes. Furthermore, the genes selected for analysis play 
important roles in the regulation of the overall gene ex- 
pression of the virus, including the iel80, the epO (and 
their antisense transcripts, the LAT and the AST, re- 
spectively), the virion host shut-off {vhs), and the ul54 
genes. For each viral gene, a minimum of 3 parallel 
replicates were performed in order to achieve statistical 
reliability. Immortalized PK-15 cells were infected with 
either the wild-type {wt) or usl gene-deleted (usl-KO) 
PRV, using a high multiplicity of infection [MOI; 10 
plaque-forming units (pfu)/cell]. The low-dose infections 
produce a much finer resolution of the cycle threshold 



(Ct) values for the transcripts than in the case of high- 
dose infections; however, in the former case a large pro- 
portion of the cells remain uninfected, which allows the 
initiation of an additional infection cycle after 6 h post- 
infection (pi) by the newly generated virions [9], which 
would confuse the interpretation of the expression data 
obtained. The transcription of PRV genes was monitored 
at 9 different time points: 0.5, 1, 2, 4, 6, 8, 12, 18 and 
24 h pi (multi-timepoint analysis). Strand-specific pri- 
mers were used for the reverse transcription reaction so 
as to exclude the distorting effects of the potential cis- 
antisense transcripts that might be produced from the 
complementary DNA strands, which cannot be avoided 
on the use of other methods, such as oligo-dT- or ran- 
dom priming-based reverse transcription [25]. On the 
other hand, we found that the specificity of strand- 
specific primer-based RT is much higher than that of 
other methods [9]. In our work, we applied a modified 
version of the mathematical model described by Soong 
and colleagues for the relative quantification [26]. Specif- 
ically, we used the average of the 6-h ECt-sample values 
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Figure 2 Heatmap visualization of the ratios of the transcripts of the usl-KO and the wt viruses (Rr = RusiKo/Rwt)- The red color (Rr < 1) 
indicates a transcript level that is lower in the mutant than in the wt virus at a certain period of infection, whereas the black (Rr > 1) indicates the 
opposite: the mutant virus produces a higher amount of mRNA than the wt virus. 
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Figure 3 The average relative expression ratios {R). (A) The R values of the total genes of the mutant and wt viruses plotted against time. 
The average level of the transcripts of the mutant virus is lower in the early stage, but higher in the late stage of infection. (B) The R values of the 
different kinetic classes of mutant and wt viruses. In the mutant virus, the L genes are repressed in the early stage of infection, whereas the E 
gene expressions are enhanced in the late stage of infection. 



of the wt PRV for each gene in both the wt and the mu- 
tant backgrounds, as controls, which were normalized to 
the average of the corresponding 28S values (ECt- 
reference). The 28S RNA gene was used as a reference 
gene since the ribosomal RNAs are not substrates of 
VHS ribonuclease [27]. We used a selected Ct value (at 
6 h in our system) as control for the comparability of 
the relative copy numbers of a transcript at different 
time points. The relative amounts of the transcripts of 
different genes cannot be compared directly due to the 
variation in the primer efficiencies in both RT and PGR. 
However, the use of multi-timepoint qPCR analysis 
allows a comparison of the transcription kinetics of 



different viral genes throughout the whole period of in- 
fection. Furthermore, this method allows a comparison 
of the same genes in the two genetic backgrounds (Rr = 
RusiKO),/Rwt)> and of the mRNAs and the complementary 
antisense transcripts in the case of the epO/LAT and 
iel 80 /AST pairs. A high R^. value indicates an inhibitory 
effect of the ICP22 protein on the transcript level of a 
particular gene in the wt virus. Conversely, a low R^. 
value indicates a stimulatory effect on the gene expres- 
sion. We applied the same logic for the interpretation 
of the data obtained as is used in other knockout 
organisms, i.e. the normal role of the usl gene is 
considered to be the opposite of that of the phenotype 
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Figure 4 The average Rr values {Ry=R usI-ko/R wt)- The thick 
solid line depicts the % values of the total genes, and the thinner 
lines the % values of the different kinetic classes of viral genes. 

V J 

caused by the mutation of this gene. In other words, an 
elevated expression of a gene in the usl-miW mutant is 
indicative of a suppressive effect of the ICP22 protein on 
the expression of this gene. 

Confirmation of the mutant genotype 

The mutation in the usl gene was confirmed by PGR 
amplification of the DNA sequences containing the mu- 
tation, followed by pyrosequencing. The mutation of the 
usl gene was rescued, and this was followed by growth 
analysis in order to confirm that the altered kinetic 
properties can be solely explained by the abrogation of 
the usl function (data not shown). Besides these ana- 
lyses, we also carried out more precise techniques for 
this purpose. Thus, we compared the rates of increase of 
viral DNA during the first twelve hours of viral infec- 
tion. We observed similar dynamics in the growth rates 
of the DNA of the wild-type and rescued viruses, while 
both differed significantly from those of the usl -mutant 
virus (Additional file 1). In addition, we compared the 
transcription kinetics of the two most important 
transactivator genes, the iel80 and epO genes of pseudo- 
rabies virus. This revealed that the kinetics of the 
rescued virus resembled that of the wild-type virus, but 
differed significantly from that of the mutant virus 
(Additional file 2 and Additional file 3). 

Transcription of the us 7 gene in the wild-type PRV 

We compared the expression kinetics of the usl gene of 
the wt virus under the low- and high-dose infection 
conditions (Figure 1). Infection was analyzed for 6 h in 
the low-MOI infection, and for 24 h in the high-MOI 



infection. The data revealed that in the low-MOI infec- 
tion the usl gene behaved as an L gene [9], since the 
amounts of its transcripts started to rise rapidly from 
4 h pi. This is in contrast with the situation of the high- 
dose experiment where this gene exhibited typical E 
characteristics in the early phase of infection, since its 
expression began to rise steeply from 2 h pi. Thus, there 
was a 2-h shift in the expression kinetics between the 
two experimental conditions up to 6 h pi. Another un- 
usual feature of the usl mRNA kinetics was the rapid 
drop in the transcript level after 6 h pi, which is not typ- 
ical in any of the kinetic classes of the viral genes. 

The impact of the usl gene mutation on the viral 
transcriptome 

The R values of individual genes were calculated for 
both wt and usl-KO mutant viruses, and compared by 
calculation of the ratios (R^. = RusiKo/^wt) for each time 
point (Figure 2). The average R values were calculated 
and plotted for the total genes (R) (Figure 3A; Additional 
file 4A) and for each kinetic class of genes {Re> R e/l and 
Rl) in both genetic backgrounds (Figure 3B; Additional 
file 4B). Additionally, the mutant and wt viruses were 
compared by calculating the average R^. values for the 
total genes ( Rr ) and separately for each kinetic class 
{R r-E^Rr-E/L ^nd R r-i) (Figure 4B). Figures 3A and 4 
reveal a significant decrease in the amounts of 
transcripts in each kinetic class of genes in the first hour 
of infection (Figure 2). The average levels of transcripts 
were not decreased significantly or at all in a few genes 
(e.g. ul21, ul24, ul53 and ull9), but in most genes, and 
especially the E genes, the fold of the decrease was close 
to 2 orders of magnitude at this very early stage of infec- 
tion. However, the average level of the E transcripts of 
the mutant virus reached approximately the same level 
as that of the wt virus in the interval 1-6 h pi. This was 
in contrast with the L genes, which on average were still 
expressed at a significantly lower level during the first 
8 h of infection. The average amounts of the E/L 
transcripts falls between those of the E and L genes. The 
iel80 gene was also expressed at a lower level in the first 
6 h of infection in the mutant virus. The deletion of the 
usl gene led to a significant reduction in the average ex- 
pression of the L genes and to a lesser reduction in the 
expression of the E/L genes in the first 6 h pi, but the E 
genes were expressed at a slightly higher level in the mu- 
tant than in the wt background in the period 2-6 h pi 
(Additional file 4B). The average rate of transcription of 
the IE, E and E/L genes from the mutant viral genome 
exceeded the rate of transcription from the wt genome 
after 6 h pi (Figure 4). This was in contrast with the 
transcript levels of the average L genes of the mutant 
virus, which remained below or equal to that of the wt 
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virus, except at 18 h pi. These data suggest that the 
ICP22 protein exerts an inhibitory effect especially on 
the transcription of the iel80 and E genes following the 
initiation of viral DNA synthesis. The usl gene deletion 
exerted the highest impact on the expression of the epO 
gene, an important regulator of PRV gene expression 
[28]: the amount of epO transcripts was ~ 10 times 
higher in the cells infected with the mutant than that for 
the wt virus within the period 8-18 h pi. This result was 
the opposite of that described in the HSV, where the 
level of ICPO mRNA was significantly lower in the KO 
virus [16]. The expressions of the us3, ul22, ul43, and 
ulS2 genes were also significantly elevated at certain 
time points in the late stage of infection (Figure 2). The 
ulSy ulSl and LAT genes were the only examined PRV 
genes whose expressions were always lower in the usl- 
KO virus than in the wt virus. Overall, the above data 
suggest that the ICP22 protein exerts a selective effect 
on the expressions of PRV genes belonging in different 



kinetic classes. ICP22 (possibly of tegument origin) 
appeared to have a significant stimulatory effect on the 
general gene expression of the virus at both 30 min and 
1 h pi [PRV ICP22 has been shown to be localized in the 
viral tegument layer [29]. Though to a lower extent, its 
stimulatory effect continued up to 6 h in the E/L and L 
genes, but the effect became neutral or slightly inhibi- 
tory in the later stages of infection in these kinetic 
classes of genes. The expressions of the E genes 
exhibited different profiles on average in the mutant 
background: no discernible effect of the mutation within 
the 1-6 h pi period, followed by a profound inhibitory 
effect in the later phases of infection. Our results may 
resolve the debate as to whether ICP22 represses the E 
genes [14] or enhances the L genes [15]. The data dem- 
onstrate that the rates of transcription of the L genes in 
the mutant virus are repressed (and therefore they are 
selectively enhanced in the wt virus) in the very early 
phase of infection, whereas the rates of transcription of 
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Figure 6 Synthesis of viral DNAs. The replication of viral DNA was 
monitored in the wt and the usl-KO viruses through the use of real- 
time RT-PCR. A 4-h delay was observed in the onset of 
DNA replication. 



the E genes are selectively enhanced (and therefore they 
are repressed in the wt virus) in the late phase of 
infection. 



Comparison of the effects on gene expression in three 
mutant PRV strains 

We compared the gene expression profiles {R r values) of 
three mutant viruses: the vhs-KO [30], epO-KO [28] and 
usl-KO viruses (Figure 5). On a broad scale, the usl and 
vhs gene deletions appeared to produce a similar overall 
expression pattern (though the usl deletion exerted a 
more marked effect), whereas epO deletion led to an over- 
all expression profile complementary to those of the usl- 
KO and epO-KO viruses. Figure 5 reveals that the usl and 
vhs gene deletions resulted in down-regulation of the PRV 
genes in the E phase and in up-regulation of these genes 
in the L phase, whereas a reverse expression pattern was 
observed in the epO mutant background. Furthermore, in 
both the usl-KO and the vhs-KO viruses, the most 
affected gene was the epO^ which was significantly up- 
regulated from 2 h pi. These results suggest a fundamental 
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Figure 7 Heatmap visualization of the ratios of the transcripts of the usl-KO and the wt viruses normalized to the relative copy 
number of the DNAs (Rm = Rusiko^Rdna-usiko/Rwi^Rdna-wi)- The values were very low in the first 2 h pi. There was a drastic change in the 
period 4-6 h pi, when the values became high. In the last stage of infection, the gene expressions became the same in the two viruses. 
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Figure 8 The average transcription rates are drastically decreased after the onset of DNA replication. The transcriptions from the genome 
of both genotypes are virtually shut down following the initiation of DNA synthesis. The transcription from the usl-KO genome is suppressed to a 
lower extent than that of wt virus. 



role of the EPO transactivator in the formation of the glo- 
bal gene expression profile of PRV. 

The expression of the ie180 gene is correlated with the 
expressions of the PRV genes in the usl-KO virus 

We previously reported that the expression of the iel80 
gene, the major transactivator of the wt virus, was 
uncorrected with the expressions of the remaining genes 
[30]. When we investigated this relationship in the usl-KO 
virus by using Pearson correlation analysis, very high 
correlations emerged between the iel80 transcripts and 
the transcripts produced by all three kinetic classes of the 
PRV genes (Additional file 5A). The correlation was 



especially high between the iel80 and the E transcripts. 
Interestingly, we observed a similar effect in the vhs-nuW 
mutant as concerns of the correlation between the iel80 
and the other viral genes [30]. 

The initiation of DNA replication is delayed in the usl-KO 
as compared with the wt virus 

Our real-time RT-PCR investigation of the kinetics of 
DNA replication of the two viral strains demonstrated 
that the amount of viral DNA was slightly less at 1 h pi 
than at 30 min pi in both viruses (Figure 6). We can only 
speculate about this phenomenon. It is possible that the 
cellular DNases digest some proportion of the infecting 
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viral DNA molecules. DNA synthesis was initiated be- 
tween 2 and 4 h pi in the wt virus. There was a decrease 
in the rate of multiplication of the DNA within the inter- 
val 4-6 h pi, which might be due to the switch between 
the theta-type and rolling circle-type of replications. A 
high rate of amplification of DNA molecules was observed 
between 6 and 12 h pi, followed by a slow decrease in the 
copy number of the wt viral DNA. The decrease in the 
amount of viral DNA in the late stage of viral infection is 
explained by the egress of mature virions from the 
infected cells. The onset of the DNA replication of the 
mutant virus exhibited a significant delay relative to that 
of the wt virus. The usl-KO DNA started replication 4 h 
later, and the rate of production of viral genomes was 
lower than that of the wt virus within the period 6-12 h 
pi. The copy numbers of the viral genomes were almost 
exactly the same at their peaks in the two viruses, but the 
maximum amount was reached 6 h later in the mutant 
virus. The mechanism whereby ICP22 affects the DNA 
replication remains to be determined, but it may be 
associated with the effect of the usl mutation on the delay 
of the onset of the general gene expression. Since the 
expressions of several key viral factors, such as the IE180 
and the EPO, were profoundly changed in consequence of 
the mutation, the effect of ICP22 on the DNA replication, 
at least in part, must be indirect. 

Investigation of the expression kinetics of PRV genes in 
the mutant and wt viruses through the use of R values 
normalized to the relative copy numbers of the DNA 

We analyzed the viral gene expressions by normalizing the 
amounts of transcripts to the relative copy number of the 
DNA molecules, which allows a comparison of the tran- 
scriptional activity with the corresponding copy number 
of the DNA at different time points. The ratios of the 
normalized R values (Rm = Rn-usiKo/Rn-wt) indicated that 
the level of transcripts in most genes was lower in the mu- 
tant than in the wt virus within the period 0.5-1 h pi with 
two exceptions {ul21 and ul24 at 30 min pi). This ten- 
dency was substantially reversed at 4 h pi, when the genes 
of the mutant virus were more active than those of the wt 
virus (Figure 7). Some early genes {epO, ul9, ul30, and 
ul52) produced > 100 times more transcripts from a single 
copy of the genes in the mutant than in the wt virus. 
There was not such a large difference between the non- 
normalized genes of the two viruses. The switch from low 
to high amounts of transcripts occurred much later in the 
mutant virus, the exact time of the switch varying from 
gene to gene (Figure 2). The average expression from a 
single gene was very low at 30 min and 1 h pi in the usl- 
KO virus as compared with the wt virus, but it had be- 
come > 50 times higher (RusiKo/Rwt) by 4 h pi (Figure 8A). 
This value was especially high for the average E genes 
(101.5), lower for the E/L genes (53), and much lower for 



the L genes (15.2). This phenomenon obviously resulted 
from the delayed onset of DNA synthesis in the usl-KO 
virus, which suggests that DNA replication exerts a severe 
constraint on the transcription. Indeed, the gene transcrip- 
tion rate was profoundly decreased after the initiation of 
DNA synthesis (Figures 8 A and SB). In the wt virus, the 
average rate of gene transcription was 33 fold less at 4 h pi 
than at 2 h pi (R4h/R2h)> and remained low or even became 
lower. In the mutant virus, the decrease in the rate of tran- 
scription started later (from 4 h pi), and occurred grad- 
ually, not suddenly as in the wt virus. The steepest drop in 
gene expression in the interval 4-6 h pi did not coincide 
with the highest rate of DNA synthesis, at 6 to 8 h pi, 
which suggests that the inhibition of overall gene expres- 
sion is not directly controlled by DNA replication itself, as 
it was suggested by Huvet and colleagues [31]. From data 
obtained in the analysis of human cells, those authors 
proposed novel gene expression regulatory mechanism, 
based on the collision of the transcriptional and DNA rep- 
lication machineries. However, those results could not be 
reproduced by Necsulea and colleagues [32]. Although, it 
is generally conceived that the switch in expression from 
the E to the L kinetic class of genes is triggered by DNA 
replication. Figure SB clearly shows that the maximum ex- 
pression as regards the individual L genes occurs at 2 h pi, 
i.e. before the onset of DNA replication, which casts se- 
vere doubt on the validity of this idea. 

Normalization of the R values to the relative DNA copy 
number leads to the ielSO gene correlating with the 
other viral genes in the wt background 

Normalization of the R values of the cDNAs to the copy 
number of the viral DNAs leads to a very strong correl- 
ation between the iel80 and the other genes in both gen- 
etic backgrounds (Additional file 5B). It is noteworthy that 
the strength of the correlation between the iel80 and the 
different kinetic classes of genes is the opposite that in the 
non-normalized case. These data suggest that the 
synchronism in the transcription between the IE 180 
transactivator and the other genes is indicative of a real 
correlation. However, since the DNA replication imposes 
a severe constraint on the transcription of each gene, in- 
cluding the iel80 gene, the high correlation between the 
normalized data could possibly be merely a statistical curi- 
osity without any functional significance. Additional file 6 
shows the R values of iel80, the average total viral genes 
and the average of various kinetic classes of genes 
normalized to the DNA copy number at different time 
points of infection. 

Conclusions 

Our kinetic data show that the abrogation of the usl 
function leads to a significant reduction of transcription 
in every kinetic class of genes in the first hour pi. In the 
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period 1 to 6 h pi, the L genes are selectively down- 
regulated, and the E genes are later up-regulated in the 
usl'KO background relative to the wt virus. The 
questions as to whether the ICP22 protein exerts a dir- 
ect effect upon the gene expression and, if so, at what 
levels remains to be answered. It is noteworthy that de- 
letion of the vhs gene resulted in a similar overall ex- 
pression pattern as that observed in the usl-KO virus, 
and in an expression profile complementary to those in 
the epO'KO virus. Interestingly, epO was the gene that 
was most affected in both the usl and vhs knockout 
viruses (highly elevated expressions), which might imply 
that EPO could be a common link in the determination 
of the overall gene expression profile. We observed a 
strong inhibition of transcription during DNA replica- 
tion. The question arises of whether the process of DNA 
replication itself could exert this inhibitory effect. In the 
mutant virus, the main inhibition in gene expression 
occurs between 4 and 6 h pi, while DNA synthesis 
exhibits the highest rate 2 h later. This suggests that the 
decrease in gene expression is not, or not only, a result 
of the potential collision of the transcription and DNA 
replication machineries. We have earlier shown that dis- 
ruption of the function of the vhs gene of the virus 
results in synchronization of the gene expression profile 
of the iel80 and the remaining PRV genes [30]. We 
obtained very similar results with the usl gene-deleted 
PRV: the expression of the iel80 gene correlated with 
those of the other PRV genes in the mutant virus, which 
was not the case for the wt virus. The question may be 
posed of whether this correlation is directly associated 
with the disruption of the usl gene function, or is rather 
caused by the delay of DNA replication. 

Methods 

Cells and viruses 

Monolayer cultures of porcine kidney epithelial cells 
(PK-15) were used for the propagation of the pseudo- 
rabies virus. Cells were grown in DMEM (Sigma 
Aldrich), supplemented with 5% fetal bovine serum 
(Gibco) and 80 [ig of gentamicin per ml (Invitrogen) at 
37°C in the presence of 5% CO2. The Kaplan wt strain of 
PRV was used as the parental strain for the generation 
of the usl -null mutant virus (usl-KO), 

Construction of the usl gene knockout virus 

The Ka-usl'KO {usl-KO, in short) virus was generated 
as follows. As a first step the BamHl-10 fragment of 
PRV was isolated from the gel, and was then subcloned 
to pRL525 [33], resulting in the generation of pRL525- 
BIO. This plasmid was used as a template for the PGR 
amplification of the two arms of the flanking sequences 
providing homology with the target viral genomic re- 
gion. A green-fluorescent protein (GPP) gene expression 



cassette (Glontech) was inserted into the unique Ed 13611 
site of the targeting sequence, resulting in pUSl-gfp, 
which was used as the transfer plasmid for the gener- 
ation of the knockout virus. The linearized transfer plas- 
mid was transfected along with the purified wt viral 
DNA into PK-15 cells. The recombinant virus was 
generated by homologous recombination, then isolated 
and plaque-purified on the basis of the fluorescence. 
Rescued viruses were generated by using pUSl as a 
transfer plasmid, which was co-transfected with the 
purified DNA of usl-KO to PK-15 cells. The revertant 
viruses were selected on the basis of the non-fluorescent 
plaque phenotype. Both mutant and rescued viruses 
were checked with DNA sequencing. 

Infection 

The virus stock used for the experiments was prepared 
by infecting PK-15 cells with low-dose viruses, followed 
by incubation of the cells until a complete cytopathic ef- 
fect was observed. To assess the effect of the usl gene 
deletion on the transcription kinetics of PRV, rapidly- 
growing semi-confluent PK-15 cells were infected with 
the wt or usl-KO virus at a high multiplicity of infection 
[MOI; 10 plaque forming units (pfu)/cell], and incubated 
for 1 h, which was followed by removal of the virus sus- 
pension and washing of the cells with PBS. Subse- 
quently, fresh culture medium was added to the cells, 
which were further cultivated for an additional 0.5, 1, 2, 
4, 6, 8, 12, 18 or 24 h. 

Isolation of RNAs 

Infected or non-infected PK-15 cells were washed with 
PBS and harvested for RNA purification. Total RNA was 
isolated by using the NucleoSpin RNA II Kit (Macherey- 
Nagel GmbH and Co. KG) as recommended by the sup- 
plier. Briefly, harvested cells were collected by low-speed 
centrifugation, and lysed in the buffer included in the 
kit. Samples were treated with RNase-free rDNase solu- 
tion (included in the Kit) to remove potential genomic 
DNA contamination. As the next step, the potential re- 
sidual DNA contamination was removed by using Turbo 
DNase (Ambion Inc.). Subsequently, RNA samples were 
eluted in RNase-free water (supplied with the kit), 
resulting in a total volume of 60 [A of RNA solution. 
RNA concentrations were measured spectrophotomet- 
rically in a BioPhotometer Plus (Eppendorf). The RNA 
solution was stored at -80°C until use. 

Reverse transcription 

Total RNA samples were reverse-transcribed by using 
gene-specific primers, and Superscript III reverse tran- 
scriptase (Invitrogen) as described in our earlier reports 
[9,34]. Briefly, RT mixtures containing total RNA, pri- 
mer. Superscript III enzyme, buffer, dNTP mix and 
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RNase inhibitor (RNAsin, Promega) were incubated at 
55°C for 1 h. The amplification of the first-strand cDNA 
synthesis was terminated by keeping the samples at 70°C 
for 15 min. The cDNAs were diluted 10-fold with 
nuclease-free water (Promega Corp.) and the solutions 
were stored at -80°C until use. 

Real-time PGR 

SYBR Green-based (Absolute QPCR SYBR Green Mix, 
Thermo Scientific) quantitative real-time PCRs were 
performed on the first-strand cDNAs in a real-time PGR 
cycler (Rotor-Gene 6000, Corbett Life Sciences), as 
described in our previous studies [9,30]. The specificity 
of the reverse transcription and the PGR reactions was 
ensured by using no-RT, no-primer, and no-template 
controls. The accuracy of sampling was guaranteed by 
using 28S rRNA of swine as loading control. The specifi- 
city of the PGR products was confirmed by melting 
point analysis, PAGE and/or DNA sequencing. 

DNA sequencing 

We used pyrosequencing with a Pyromark Q24 pyro- 
sequencer (Qiagen) to validate the mutation of the usl 
gene and the specificity of the PGR products obtained in 
the kinetic experiments in the event of doubt. 

Data analysis 

The relative expression ratios (R values) of the cDNAs 
of the PRV genes were calculated via the following for- 
mula: 

\psample6h ) _ y-'ref 6h ) 

[^sample ) \ ref ) 

The cDNAs were all normalized to the cDNAs of the 
28S rRNAs of swine [9] by using the Gomparative 
Quantitation module of the Rotor-Gene 6000 software 
(Version 1.7.28, Gorbett Research), which automatic- 
ally sets the thresholds and calculates the efficiency of 
PGR reactions. We used the average 6-h EGt values of 
the "samples", with those of the "references" as 
controls, as in our earUer works [9,30,34]. The R values 
of the viral DNAs were calculated similarly; the 6-h 
EGt values were taken as the control, and porcine 28S 
rRNA gene was used as the reference. The effect of de- 
letion of the usl gene on the global gene expression 
was calculated by using R^ the ratio of the R values of 
the usl mutant and the wt viruses (R^. = RusiKo/Rwt)> 
where Rusiko and R^t are the R values of a particular 
gene at a given time point in the usl -KG and wt gen- 
etic background, respectively. All data were analyzed 
by using the average and the standard deviance 
functions of Microsoft Excel. Pearson's correlation 
coefficient was calculated for the analysis of the 



correlation between the expression kinetics of the 
genes, using the following formula: 

n 

Y^{Xi-X){Yi-Y) 

r = ^ {n-i)s s [^^]* normalized R 

values were calculated by dividing the appropriate R 
value of the cDNA by the R value of the viral DNA 
measured at the same time for the same sample. The 
Pearson correlation coefficient is a number between -1 
and +1 that measures the linear relationship between 
two variables, denoted here as X and Y, which are the R 
values of two different genes or the average R values of 
genes belonging in the same kinetic class in the same 
time interval. X and Y are the average values, n is the 
sample number, and SX and SY are the standard errors 
of the mean values for X and Y, respectively. A positive 
value for the correlation indicates a positive association, 
while a negative value indicates an inverse association. 

Additional files 



Additional file 1: Comparison of the rates of increase of viral DNA 
during the first twelve hours of PRV infection. We observed similar 
dyr^amics in the growth rates of the DNA of the wild-type and rescued 
PRVs, while both differed significantly from those of the usl -mutant virus. 

Additional file 2: Comparison of the transcription kinetics of the 
ie180 gene in the wild-type, us7-K0 and us7-rescued PRV. This 
revealed that the kinetics of the rescued virus resembled that of the wild- 
type PRV, but differed significantly from that of the mutant virus. 

Additional file 3: Comparison of the transcription kinetics of the 
epO gene in the wild-type, us7-K0 and us7-rescued PRV. This 
revealed that the kinetics of the rescued virus resembled that of the wild- 
type PRV, but differed significantly from that of the mutant virus. 

Additional file 4: The average relative expression ratios {R). This 
table shows the R values for the total PRV genes (A) and for each kinetic 
class of viral genes (B) at different time points of infection. 

Additional file 5: Correlations between the transcription of the 
iel80 gene and other PRV genes in wt and usl-KO backgrounds A. 

The viral genes are expressed in synchronism with the iel80 gene in the 
mutant virus, whereas the expressions are not correlated in the lA/r virus. 
B. Correlation between the transcription of the iel80 gene and other PRV 
genes in the wt and usl-KO backgrounds with the use of normalized R 
values. The expression of viral genes becomes correlated with the 
expression of iel80 genes in the wt virus, too. 

Additional file 6: The normalized average relative expression ratios 
{R^ ) values This table shows the R^ values for the total PRV genes 
(A) and for each kinetic class of viral genes (B) at different time 
points of infection. 
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